SINGULARITIES OF AXIALLY SYMMETRIC VOLUME PRESERVING MEAN 

CURVATURE FLOW 
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^»T-> ", ABSTRACT. We investigate the formation of singularities under the volume preserving mean curva- 

ture flow. We study axially symmetric surfaces with Neumann boundary conditions and prove that 
under an additional lower height bound on the boundary of a specific region, the first singularity that 

CNI ■ 

forms is of type I. 

^P ■ 1. Introduction 

CN 

Let M n be an n-dimensional manifold and let xo : M n — >• R n+1 be a smooth immersion of M n 

into R n+1 . Consider a one-parameter family of smooth immersions x t : M n — > M. n+1 with M t = 

W . xt (M n ) and x t = x(«, t) satisfying 

^: d 

(1.1) -x(Z, t) = -H(l, t)v(l, t), l€M n ,t>0. 

By i/(l, t) we denote a designated outer unit normal of M t at x(7, t) (outer normal in case of com- 
pact surfaces without boundary), and by H(l, i) the mean curvature with respect to this normal. 
Here the hypersurfaces M t are evolving by mean curvature. 



If the evolving compact surfaces M t are assumed to enclose a prescribed volume V the evolution 
equation changes as follows: 



in 



CO 

(1.2) jx(l,t) = -(H(l,t)-h(t))v(l,t), lGM n ,t>0, 

where h(t) is the average of the mean curvature, 



h(t) 



hl t d 9t 

and gt denotes the metric on Mt. This flow is known to decrease the surface area while the enclosed 
volume remains constant. 

Huisken [15] proves that uniformly convex, compact surfaces become asymptotically spherical un- 
der mean curvature flow. Grayson [14] proves that mean curvature flow makes smooth embedded 
curves in the plane shrink to a point, becoming spherical in the limit. Ecker and Huisken [10] prove 
that entire graphs over W 1 "flatten out" with time. In [1] Altschuler, Angenent and Giga study 
mean curvature flow of surfaces which are axially symmetric and prove that immediately after a 
singularity the surface becomes smooth and thus continue to study the flow after singularities. In 
[19] Huisken and Sinestrari avoid the formation of the singularity by performing surgery to the 
surface before the singularity develops. In this procedure they remove the part of the surface which 
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becomes singular, and then patch up the remaining surface by a smooth cap, and afterwards let the 
surface evolve once again. 

If a uniformly convex, compact object without boundary flows by mean curvature while keeping 
the enclosed volume constant, Huisken [16] proves that it converges to a sphere. The first author 
[4] proves that an axially symmetric hypersurface between two parallel planes, which encloses a 
sufficiently large volume, having Neumann boundary data, converges to a cylinder under volume 
preserving mean curvature flow. In [5] she proves that type I singularities of axially symmetric 
volume preserving mean curvature flow are self-similar and asymptotically cylindrical. Escher and 
Simonett [12] investigate the volume flow near spheres. They prove that there exist global solutions 
to this flow starting from non-convex initial hypersurfaces. Cabezas-Rivas and Miquel [7] study the 
volume flow in Hyperbolic space. 

As the initial surface we choose a compact n-dimensional hypersurface Mq, with boundary dM$ / 
0. We assume Mq to be smoothly embedded in the domain between the two parallel planes x\ = a 
and x\ = b , a, b > 0. Here, except for the volume constraint, we have a free boundary. We con- 
sider an axially symmetric hypersurface contained in the region between the planes x\ = a and 
x\ = b. Motivated by the fact that the stationary solution to the associated Euler Lagrange equation 
of an energy minimizing liquid bridge satisfies a Neumann boundary condition, we also assume 
the surface to meet the planes at right angles along its boundary. We study the first singularity that 
develops under this flow. 

We assume a lower height bound detailed in Assumption 7.4 to prove the following theorem. 

Theorem 1.1. Under the Assumption 7.4, when the first singularity develops the second fundamen- 
tal form \A\ 2 satisfies 



max.\A\ z < 



2 - C 



Aft ' ' - T - t ' 
for all t <T, where T is the blow up time and C , is a constant. 

In [17] Huisken obtains type I singularities in the mean curvature flow setting for axially symmetric 
hypersurfaces with positive mean curvature. In our case, we have no conditions on the curvature, 
but we have a lower height bound on the boundary of a specific domain. 

The paper is organized as follows: 

In Section 2 we introduce some notation. 

In Section 3 we discuss parabolic maximum principles for non-cylindrical domains. As we use the 
parabolic maximum principles for specific regions of the hypersurface, our domains change with 
time, and therefore may not be cylindrical. When applying the maximum principles to a specific 
region of the hypersurface, it does not help to consider a periodic hypersurface by reflecting it along 
the planes x\ = a and x\ = b, which alleviates the need to take the boundary data into account. 
Hence the boundary data of the specific region plays an important role in our case. This version of 
the parabolic maximum principle is particularly useful as it rules out some specific regions of the 
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parabolic boundary from attaining a maximum. 

In Section 4 we compute the evolution equations and introduce different regions used in the rest of 
the paper. These different regions give us a better understanding of the hypersurface and later we 
prove that the singularity can only develop in certain regions. 

In Section 5 we compute the height, gradient and curvature estimates. In particular we prove that a 
specific region has a lower height bound for all time t < T and also the second fundamental form 
| -A | is bounded in that region. We also prove that the mean curvature is bounded from below on the 
entire hypersurface. 

In Section 6 we discuss the linear and non-linear versions of the Sturmian Theorem and its appli- 
cations in volume flow. We study the behaviour of zeros of particular functions and discuss the 
implications of the zeros. 

In Section 7 we prove Theorem 1.1 by studying the different cases in which a singularity can 
develop. 

Acknowledgment. The second author thanks Prof. Gerhard Huisken and Dr Ben Andrews for helpful discussions; and 
for supporting visits to the Max Planck Institute for Gravitational Physics and ANU respectively. 



2. NOTATION 

We follow Huisken's [17] and Athanassenas' [4] notation in describing the n-dimensional axially 
symmetric hypersurface. Let po : [«> b] — > R be a smooth, positive function on the bounded interval 
[a, b] with p' (a) = p' {b) = 0. Consider the n-dimensional hypersurface Mq in R n+1 generated by 
rotating the graph of po about the xi-axis. We evolve Mq along its mean curvature vector keeping 
its enclosed volume constant subject to Neumann boundary conditions at x\ = a and x\ = b. 
Equivalently, we could consider the evolution of a periodic surface defined along the whole x\ axis. 
By definition the evolution preserves axial symmetry. The position vector x of the hypersurface 
satisfies the evolution equation 

—x = -(H - h)v = H + hv , 
at 

where H is the mean curvature vector. Since Ax = H, where A denotes the Laplacian on the 

surface, we obtain, 

(2.1) — x = Ax + /iz/. 

at 

Letii, . . . , i n +i be the standard basis in R n+1 , corresponding to x\, . . . x n+ \ axes andri(t), . . . ,T n (t) 

be a local orthonormal frame on M t such that 

(r/(t),ii)=0, for l = 2,...,n, and <n(t),ii) > 0. 

Let tjj = -m- € R n+1 denote the unit outward normal to the cylinder intersecting M t at the point 
x(7, t) , where x = x — (x, ii) ii. Let 
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y = (x, w) and v = (u, v) 

Here y is the height function. We call v the gradient function. We note that v is a geometric 
quantity, related to the inclination angle; in particular v corresponds to \J\ + p' 2 in the axially 
symmetric setting. The quantity v has facilitated results such as gradient estimates in graphical 
situations ( see for example [8, 10] ). 
We introduce the quantities ( see also [17] ) 

(2-2) P=(Ti,h)y~ 1 , 9= (^ii)y -1 , 

so that 

(2-3) p 2 + q 2 = y' 2 . 

The second fundamental form has n — 1 eigenvalues equal to p = — , 1 and one eigenvalue 
equal to 

fc= ^ 1,/ ' Tl > = (1+^2)3/2- 

We note that p(x±,t) is the radius function such that p : [a,b] x [0, T) — > R, whereas y(l,t) is 
the height function and y : M n x [0, T) — > R. These are two different interpretations of the same 
physical object. 

There are cases where pinch off develops in the axially symmetric setting (see [5]). We assume that 
a singularity develops for the first time at t = T < oo. 

3. Maximum Principles 

We are interested in maximum principles in non-cylindrical domains. This section is an extension 
of Ecker's ([9] Proposition 3.1) and Lumer's [21] version of maximum principles to our setting. In 
[21] the maximum principles are proved in an operator theoretic setting, which has been adapted to 
the manifold setting here. 

Let Q = M n . Let V C £1 x (0, T) be an open non-cylindrical domain. Let Q t = Q, x {t} , and 
for t / let Vt = fit H V. Let V denote the closure of V. Let Vo = $\) l~l V . The boundary of 
V is dV = V\V . The parabolic boundary Ty = dV\Qr ■ To describe the horizontal parts of the 
boundary of V in the spacetime diagram, we define the following: let Z t be the largest subset of 
Qt H dV , which is open in dV and that can be reached from "below" in V, where t is the vertical 
axis. Let Zy = Uo<«T ^* ■ Let &v = ^v\Zy ■ 

Proposition 3.1. ( Maximum Principle ) Let (M t ) t€ / T -> be a solution of volume preserving mean 

curvature flow consisting of hypersurfaces Mj = x$(fi) where Xj = x(-,t) : 17 x [0, T) — > R n+1 
and O is compact. Suppose f G (C 2 ' l (V) n C(V)) a«<i / : V — > R satisfies an inequality of the 
form 

±-A)f<{a,Vf), 



dt 

where the Laplacian A and the gradient V are computed on the manifold Mf (for the vector field 
a : V — > R n+1 we only require that it is continuous in a neighbourhood of all maximum points of 
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Figure 1 . The non-cylindrical domain V, with 5y indicated by a darker line for 

t < T. 



f ). Then 

sup / < sup / . 
v r v 

for all t e [0,T). 

Assuming f to have a positive supremum in V then 

sup / < sup / . 

V 8 V 



for all te [0,T). 



Proof. Part A. We show that 
(3.1) 



sup/(M) < sup/(/,t) 

V T v 



Let / = / - eit , where ei > . It holds /(/, 0) = f(l, 0) and on F v , f(l, t) < f(l, t). We note 
that 



f = |- £1 , A/ = A/. V/-V/. 



Therefore 



d_ 



A-a-V /<0 



At an interior point of V, where for the first time maxy t / reaches a value larger than sup r /, the 
standard derivative criteria at the local maximum say 

d[ df__ d 2 f 

dt ~ dx; dxidxA ~~ 



J% u^ 3 



As 
(3.2) 



M = F \ -M- ~ r?# I and V/ = g «*L*L 



CLJL'liliJb j (Jjjbfe 



diJb j (XJb £ 



by choosing the parametrization such that gij = 5ij at the point that corresponds to the interior 
maximum we have 

'|-A-„.v)/>0. 

This is a contradiction. Hence f(l,t) is bounded by the values of sup r / at all times. Therefore 

sup (f(l, t) - eit) < sup / < sup / , 

v r> r v 



6 MARIA ATHANASSENAS AND SEVVANDI KANDANAARACHCHI 

sup /(/, t) - e-yT < sup (f(l, £) - eii) < sup / , 

V V r v 

sup /(/,£) < supf + eiT for all ei>0, 

v r> 

giving us 

sup /(/,£) < sup/, 

v ry 

which completes Part A . 

Part B. For this part we suppose that / has a positive maximum in V. We will prove by contradic- 
tion that 



(3.3) 



sup / < sup / . 

V 5 V 



Suppose that (3.3) does not hold. As supy /(/, £) < sup r ^ / this can only happen if the maximum 
of / is achieved on Zy. In particular the maximum of / can only be achieved at an interior point 
of Zy for (3.3) to be contradicted. We denote by ^ max the union of Z t 's on which the maximum is 
achieved. Let £* denote the first time that the maximum is achieved on Zy. Let 



K = {(l,t)EV:f(l,t) 



sup /} . 

v 



t = T 



t = 




Figure 2. If the maximum is achieved on Zy 

(We note that K n dV is non-empty as the maximum is achieved in Z n 
Zy. ) Therefore there exists a (3 > such that 



and also that KndV C 



sup / > j3 > sup 

V ' T v \Z m , 



/• 



As Sy C (rV\Z max ) we have 



Define 



sup 

ry\z mi 



> sup / . 

<v 



K fi = {(l,t)€V:f(l,t)>p}^<b. 

Here Kp may not be a connected set. It holds that Z tt C Kp. We work with the connected 
component of Kp, which has Z tt as a part of its boundary. Let < €2 < supy / — /?. As V is open 
there exists (^2^2) £ V (depending on €2) such that f(h,t2) > supy / — €2- We take £2 to be the 
earliest time such that /(/, t) > supy / — £2 is satisfied. Take £3 G (£2, £*) and choose a smooth 
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function <j) : [0, T] ->• R , such that < <j) < 1 , <j)' < on (t 2 ,t 3 ) , and = 1 on [0, t 2 ], and <f> = 
on [£ 3 ,T]. The set 

is a compact subset in V. The set K'o is in the interior of V, thus K^ n Zy = . As (<f)f)(l 2 , t 2 ) = 
f(h,t 2 ), and as the maximum of / is achieved in Z tt , and as <fi > , the supremum of 4>f must be 
at least as big as f(l 2 , £2)- Also 



max </> f = max < 



as = on [£ 3 , T]. Hence 



max 4>f> f{l 2 ,t 2 ) >sup/-e 2 >/3> 0. 



On the other hand, V\K' p = (V\Kp) U \Kp\K'A . On {V\K p ) n {(/, i) G F : / > 0} we have 
0/ < / < P, and on (V\Kp) n {(/, t) G F : / < 0} we have 0/ < /? . On i^V^ , <j> = 0. As a 
result 



(3.4) sup <pf < (3 < maxc 



v\n k 'p 



Therefore if we denote by (l^,,t^) the point at which (</)f)(l^,t^) = sup v (4>f) , we can see that 



-|(0/) = <Af + <A7, A(0/) = <^A/, V(<f>f) = tj>Vf, 



we have 



As </>'(/</,, t^) < and /(/«£, £</,) > we obtain 



|-A-cv) W )<0. 
By using Part A with (f>f replacing / we have 

supO/) < swp((/>f) . 
v r v 

But this is a contradiction as supy (<f>f) = (</>/)(/</>,£</>) with (l^,t,p) $. Ty , and because Ty C 
V\l£a and (3.4) holds. Therefore our original assumption is wrong. Hence a maximum of / does 
not occur in Zy, that means K n dV <£ Zy . Therefore we conclude that (3.3) is true. 

□ 

4. Evolution Equations and Different Regions of the Hypersurface 

4.1. Evolution equations. The computation of the evolution equations is very similar to that of 
([17] Lemma 5.1 ) and ([4] Lemma 3 ). 

Lemma 4.1. We have the following evolution equations: 

(i) j§<x,*i) =A(yL,i 1 ) + hqy; 
(ii) f t y = Ay - 2=i + hpy ; 
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(iii) jL q = Aq + \A\ 2 q + q((n - l)p 2 + (n - 3)<? 2 - 2fcp) - frpg ; 

(iv) j^p = Ap + |,4| 2 p + 2g 2 (/c - p) - Zip 2 ; 

(v) f t k = Ak + \A\ 2 k - 2(n - 1)^ 2 (A; - p) - ZiA; 2 ; 

(vi) |iT = Air+(JJ-fc)|A| 2 ; 

(vii) f t \A\ 2 = A\A\ 2 - 2\VA\ 2 + 2\A\ 4 - 2/rif ; 

(viii) iv = Av- \A\ 2 v + (n - 1)^ - §|V^| 2 ; 

where ^ = g v g kl g mn h ik hi m h nj . 



Proof. The first identity is immediate from (2.2) and (2.1). 

Since |x| = y , we write y as y = (|x| 2 — | (x, i\) j 2 ) 1 / 2 as . To obtain the second identity we 

compute 

(4-1) — -y = y~ 1 ((H,x) - (x,ii) (H,n) + h(u,x) - hqy (*, h)) , and 

at 

Ay = y- 1 (n + (x, H) - (x, h) (H, h) - \ ( Tl ,h) | 2 - \Vy\ 2 ) . 

Using the relation 

I fa, *i> | 2 = yV = 1 - y 2 q 2 = 1 - |Vy| 2 , 
we find 

— - A J y = jT 1 (h (u, x) - hqy (x, ii) - (n - 1)) . 

Then (ii) follows after considering qy = {v,i\) and 

(z/,x) - (Mi)(x,ii) = (y, (x,ii)ii +x) - (i/,ii)(x,«i), 

(4.2) = (i/,x) = |x|(i/,w) = yt^ 1 = py 2 . 

Also using (4.2) and (4.1) we derive an alternate evolution equation for y which will be of use later. 

-r + y = 0,x> -(x,ii)(i/,ii)) , 

o(t y 

it y. 

(4.3) = py 2 = -(H- h)py . 



y 



d 



For (iii) we note that %v = VH (as shown in [15], Lemma 3.3), and using Codazzi's equation we 
obtain 



We note that 



— (u,i 1 ) = A(u,i 1 ) + \A\ 2 (u,i 1 ) 



d . d dy dq 

— (v, i\) = — (ay) = q—. — \-y— and 

dt x ' ' dV Hy ' H dt y dt 



A (i/, ii) = A{qy) = qAy + yAq + 2{Vq, Vy) . 
Combining this with (ii), we find 

y— = (qAy + yAq + 2{Vq, Vy)) - q (Ay - (n - l)y _1 + hpy) + |.4| 2 gy , which gives 
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(4.4) 

J = Aq + -(Vq, Vy> + \A\ 2 q + (n - 1)-| + hpq . 
at y y 2 

We know from ( [17] Lemma 5.1 ) 

(4.5) Viy = 5nqy, Vi {v,h) = kpy , V^g = 5ii (g 2 + fep) . 

Hence we obtain 

2 o 

-(V?, Vy) + (n - 1)-% = Q (-2fcp + (n - l)p 2 + (n - 3)g 2 ) . 

y y^ 

Using this in (4.4) gives the evolution equation for q. 
From (ii) and (iii) we can now compute 

d d i —1 -2\i/2 

Jt P= dt {y ~ q ) > 

2 P V y 3 dt dt 
= - < — g ( Ay h /ipy j - g (Ag + |^4| 2 g + q((n - l)p 2 + (n — 3)g 2 - 2/cp) - hpq) > , 

J ' 1 qn 1 [IT) 

- - { —Ay + qAq + — £ - |A| 2 g 2 - g 2 ((n - l)p 2 + (n - 3)g 2 - 2/cp) + /ipg 2 



p U° y" y" 

Ap = - ( ^Ay + gAg - |Vg| 2 + -^|Vy| 2 - |Vp| ; 

p \y y 

d -p = Ap + ^-A-|^| 2 ^V((n-l)p + (n-3)^-2 fc ) + V + LM 3 |Vz/|2+ |Vp| 2 



(it 1 ' *" py 4 y 2 p ' p p py 4 *" ' p 

Equation (iv) follows if we combine (4.5) and the following relations: 

(4.6) Vip = 5 iiq (p-k), \A\ 2 = k 2 + (n - l)p 2 , y~ 4 = p 4 + 2p 2 g 2 + q 4 . 

The evolution equations for H and |yl| 2 were derived in [16], while (v) follows from (vi), (iv), and 
the fact that H = k + (n — l)p. The evolution equation for v was derived in [4] and for p in [5]. D 

Lemma 4.2. We have V n u> = (1 — 8u) -t\ , where V is the covariant derivative on IR n+1 . 

Proof. For any unit vector rj = rj k i k we have 

71+1 W+l / n \ 

V^x = r/ fc V ifcX = J^ % J^ ( ft^ Xi + r ' fcX ' lj 

k=l J— 2 
n+1 n-f-1 n+1 



X] % X] ^^' = $Z %ifc = ^ ~ (^ ' ^) ^ > 

fe=l i=2 fc=2 



as T ik , the induced connection on M. n+1 is zero with orthonormal coordinates. And 

d 

77(1*1) =%ifc(|x|) =%-— |x| 

= 7] k uj k = (rj ,uj) , 
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as bj\ = 0. So we obtain 

V„0J = Vn ( Tn 

\|x| 

= — V r ,x- — 2?7(|x|)x 
|x| |x| z 

y 

Therefore for 77 = n we have 

(4.7) V Tl w = - (n - (n , ii) ii - (n , w) w) = . 
For 77 = 77 , Z 7^ 1 we have 

(4.8) V r; u; = - (rj - (t/ , ii) ii - {n ,oj)u) = -n . 

v y 



Corollary 4.3. The gradient function v satisfies the following alternate evolution equation 

— - A J v = -\A\ 2 v+p 2 v 3 - 2vk 2 {v 2 - 1) . 
Proof. First we compute an expression for |Vv| as follows: 

Vv = Ti {V ,w)~ Ti 

= -V 2 (( Vjf , w) + {v , Vjw}) Ti 

1 



□ 



-v 2 Uh ijg ^T k ,uj^ + /iy,0 



i - 6u) -Ti ) I Ti 

y 

as ViU = hijg^ T k and from Lemma 4.2 . 



v 2 hijg jk (r k , w) 



TV 



= -v 2 hijg jl (ri,w)rj 
= -w 2 /in5 n (-gy)ri, 
as g= (z^ , ii)^^ 1 and (i/,ii) = -(n ,u;) 



2, ,9 * 

■u kqyT\ = vk—Ti , as p = — 



(4.9) |Vv| 2 = v 2 £: 2 ^ . 

From Lemma 4. 1 we know 



^ A \ I 4,9 ^ 2,„ ,, 

— - A u = - 4 2 u + -r - -|Vv| 2 . 
at J y z v 



Substituting p = — , (4.9) and \ = v 2 — 1 , we obtain 



— - A ) v = -\A\ 2 v+p 2 v 3 - 2vk 2 (v 2 - 1) . 
at 



□ 
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4.2. Bounds on h. We state [5] Proposition 1.4 here. 

Proposition 4.4. (Athanassenas). Assume {Mt} to be a family of smooth, rotationally symmetric 
surfaces, solving (1.2) for t £ [0, T) . Then the mean value h of the mean curvature satisfies 

< C 2 < h < C3 , 

with c 2 and C3 constants depending on the initial hypersurface Mq. 

This is an important result that will be used repeatedly in our paper. 

4.3. Different regions of the volume flow surface. Depending on the situation, we are interested 
in different parts of the hypersurface. Therefore we subdivide the hypersurface in different regions 
as follows : 

4.3.1. The regions tit and tip Let 

ti t = {x(7, t) £ M t : H(l, t) < c 2 } and ti = (J ti t , 

t<T 

and ti t = {x(Z, t) G M t : H(l, t) > c 2 } and ti = (J ti t , 

t<T 

such that M t = ti t U ti t . In view of Proposition 4.4 and equation (4.3) we can see that in the region 

ti the height function is bounded from below. 

4.3.2. The regions ti' t , ti' t and ti". Let 

ti't = |x(/,t) e M t : H(l,t) < — J and ti' = (J ti' t , 



t<T 



ti 



' t = |x(Z, t) eM t : H(l, t) > — J and ti' = \J ti' t 



t<T 



such that M t = ti! t U ti' t . We also define 

ti" = {x(Z, t)eM t : H(l, t) < c 2 - S, S > 0} and ti" = \J ti'{ , 

t<T 

which will be used occasionally. 

4.3.3. The regions ti* + and ti*~ . Let 

tf+ = { X ( Z) t)€M t : k(l, t) > 0} n {x(/, t)GM t : H(l, t) > 0} , and 
£l k t - = {x(Z, t)GM t : k(l, t) < 0} n {x(/, t)€ M t : H(l, t) > 0} . 

4.3.4. The regions ti^ and ti\. For a constant cqq > 2, which we will choose later, let 

^{ x( M).M,:iM<y=^r} -drf-um, 
m= {x,m) e m, h > ^^} -rf-um, 

such that M t = f^ U ti\. 
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5. Height, Gradient and Curvature Estimates 
5.1. Height estimates. The first author proves in ( [4], 2A Remark (iii)) that the height y satisfies 

y<R, 

for some R > determined by the initial hypersurface Mo. 

We will show that the height function y has a lower bound in the region Cl'. 

Lemma 5.1. There exist constants c, d > such that inf^, y = infr-, y > c and infr-, y > d , 
where T^, and T^, denote the parabolic boundary ofCl' and Cl' (see figure 3) respectively. 




Figure 3. Space time schematic diagram for Cl and Cl' 

Proof. As -3L = — (H — h)py > in Cl' , the height increases in this region. Therefore 

inf y = inf y . 

We claim that infr 5 , y / . To prove this suppose infr~, y = , at a point x(/o, to) G Cl' , (see 
figure 3) where Cl' is the closure of Cl' and to ma y equal to T. If the height is zero at the point 
x (^o> *o) ^ Cl', then the height has to decrease near the point just before to- That means there exists 
a neighbourhood N of x(Zo, to) j suc h that N is "past" in time, t < to, and N C Cl , and -^ | N < . 
But this is not possible, since ^jfUi„t > , where Cl mt denotes the interior of Cl. Therefore there 
exists a constant c such that, on the parabolic boundary of Cl' , inf y > c > . 
When we consider M t as a periodic hypersurface, we have dCl' t = dCl' t , that is r^,\f£ = r^,\J7 . 
As inf Mo y/Owe have the desired result. □ 

Remark 5.2. The lower height bound similar to that in Cl' can be obtained in Cl" for any 5 > 0. 

5.2. A Gradient estimate. 

Lemma 5.3. There exists a constant c^ depending only on the initial hypersurface, such that vy < 
C4, independent of time. 
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Proof. We calculate from Lemma 4. 1 

d 2 

— (yv - c 3 t) = A(yv) (Vv , V(yv)) - yv\A\ 2 + h - c 3 . 

at v 

As h < cs we get by the parabolic maximum principle 

yv — c%t < max yv , 

M 

yv < max yv + C3T =: C4 . 

□ 



5.3. Curvature estimates. 

Proposition 5.4. There is a constant c\ depending only on the initial hypersurface, such that - < 
C\, independent of time. 

Proof. Similar to equation (19) of [17] we calculate from Lemma 4.1 

d fk\ A k 2/ /k\\ q 2 , TW/ N , hk , 

- - = A- + - ( Vp , V - ) + 2% (p - k) n - 1 p + k) + — (p - k) . 
at \p J p p\ \p J I p z p 

If - > 1 then (p — k) < 0. By the parabolic maximum principle we obtain 

k ( k\ 

(5.1) — < max I 1, max — I =: c\ . 

p \ Mo pj 



D 



Proposition 5.5. At points x(Z, t) of M t where H > we have — < max(ci, n — 1). 



fc + N 



Proof In a region or at any given point where H is positive, if k is positive as well (in Cl t ) we 
have by Proposition 5.4 — = - < c±. If k is negative, then 

k + (n- l)p> 0, 
-\k\ + (n-l)p>0, 

(5.2) iM<( n _i). 

□ 

Proposition 5.6. Forgiven e > 0, let St C M t ara<i S = Ut<T ^t ^ e a re gion such that y\s > e > 
anii i? |a5 t > 0/or allt <T . Then the norm of the second fundamental form \A\ is bounded in S. 

Proof. We proceed as in ([11], proof of Theorem 3.1) and ([4], Proposition 5) and calculate the 
evolution equation for the product g = \A\ 2 tp(v 2 ), where ip(r) = A _ r - , with some constants 
X, fi > to be chosen later and v = (u, lo)^ 1 . From the evolution equation of g we find the 
inequality 

j t - A) g < -2w>-2\<pv- 3 (Vv , Vg)- - A ^ |Vt;| 2 g-2hC ^v 2 )^^^v 2 ^\A\ 2 . 
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We estimate the second last term as in [4] using Young's inequality and obtain 

-2hC^{v 2 ) < 2h\A\ 3 <p(v 2 ) 

Note that by Lemma 5.3 v is bounded in S, v < — . We choose n > | and A > //maxi; 2 . As 

^ = (A-/^)^ We haVe 

2(n- 1) 2 .. il2 2(n-l)A 

y z y (A — /tii^J 

Therefore we obtain 

j t - AJg < - Un - g 2 + c 85 - 2A^~ 3 (Vt> , Vg) + ±/i V 2 • 

We note that h 4 ip~ 2 (v 2 ) and 2A(^w -3 are bounded as v is bounded in S. By relabeling the constants 
we obtain 

|-a) 5 <-c 9 ( 5 -^) -c 10 (Vv,Vg)+ Cll . 

For ^ > ^- + . ^ the right hand side of the last inequality is negative. From the parabolic 
maximum principle 

, ( c 8 f^T\ 

q < max max a , \- . / — , 

where Ts denotes the parabolic boundary of S. Therefore 

(5.3) |A|W)<max max|A|V(A max |A| V(A ^- + */^ 

\ «T 

As H\ds t > for all t < T , we have -y\ ds < max(ci, n — 1) =: c 1 for all t < T , by Proposition 
5.5. Therefore 

\A\ 2 \ dSt = (k 2 + (n - l)p 2 )\ dSt < (n - 1 + cVlafc , 
< (n - 1 + c^y- 2 ^ < (n - 1 + c^ 2 , 

for all t < T. Note that tp{v 2 ) > and is bounded from above as long as v is bounded, which holds 

for any points that are at a distance larger than e from the axis of rotation (Lemma 5.3). Therefore 

max^Sj | A| 2 (/?(t; 2 ) is bounded and this completes the proof. □ 

t<T 

Lemma 5.7. For x satisfying (1.2), the norm of the second fundamental form \A\ is bounded in the 
region ti'. 

Proof. From Lemma 5. 1 we know that in Ct' , inf y > c > . On the boundary of Cl' t , H = ^ > 
for all t < T. Therefore by Proposition 5.6 there exists a constant C such that |^4| 2 |q/ < C" < oo . 

□ 
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Lemma 5.8. There exists a constant C independent of time such that H(l,t) > —C 2 for all 

x(l,t) G M t . 

Proof By definition Cl' = \J t<T {x(l,t) G M t : H(l,t) < ^}. H can only be negative for x(l,t) G 
&'. But by the above result |-4| 2 |$-y < C < oo . As \H 2 < \A\ 2 we deduce 




H\ Mt >-V^C'=:-C 2 

□ 

Now we refine Proposition 5.6 to show that no singularities develop away from the axis of rotation. 

Proposition 5.9. For given e > 0, let St C Mt and S = \J t <T &t , such that y\s > e > ,for all 
t < T . Then the norm of the second fundamental form \A\ is bounded in the region S. 

Proof. The proof is exactly the same as in Proposition 5.6 up to (5.3). We compute the evolution 
equation for\A\ 2 tp{v 2 ) and end up with 

(5.4) \A\ 2 ip{v 2 ) < max max \A\ 2 cp(v 2 ), max \A\ 2 (p(v 2 ), 

\ So dS t 

\ t<T 

Let 

dS? = {x(/, t) G dS t : H(l, t) > 0} , and 

dS^ = {x(/,t) S dS t : H(l,t) < 0} , 
so that dS t = dSf U dSf. Continuing from (5.4) 

l^l^^fT; 2 ) < max I max \A\ cp(v ), maxlAI (f(v ), maxlAI ip(v ), \- . — 

So os+ dSf 2c 9 V c 9 

t<T t<T 

Here we look at the term max gs - |^4| 2 c/?(t; 2 ) , as the other terms are taken care of in Proposition 

t<T 

5.6. As H < < f on dS^ , dSf C Q.' t . From Lemma 5.7, \A\ 2 is bounded in Q.' . As dS^ C fy, 
\A\ 2 is bounded on dS^ for all t < T. As (p(v 2 ) > and is bounded from above at points away 
from the axis, max as - |^4| 2 (/?(v 2 ) is bounded and thus we get the desired result. □ 

t<T 

6. Sturmian Theorem and its Applications 

6.1. Introduction. Assume / : R 2 — >■ R to be a solution of 

(6-1) f t = a(z,t)f zz + b(z,t)f z + c(z,t)f, 

on Q = {{z, t) G M 2 : < z < 1, < t < T} with Dirichlet or Neumann boundary conditions. 
The number of zeros of /(•, t) is defined as the supremum of all k such that there exist < z\ < 

Z2 < ■ ■ ■ < Zk < I with 

f ( Zi ,t)f(z l+l ,t) < , i = 1,2, ... ,k - 1 . 

For t G (0, f) let 

Z t (f) = {xeR:f(x,t) = 0}, 
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be the zero set of /. In [2] Angenent proves the following theorem. 

Theorem 6.1. (Angenent [2], Theorem C) . Let f : [0, 1] x [0, T] — > R be a bounded solution 
of (6.1) which satisfies either Dirichlet, Neumann or periodic boundary conditions. Assume the 
coefficients of (6.1) to satisfy 

a > 0, a, aT x ,a u a z ,a zz G L°° , 

b,b t ,b z GL°°, 

ceL°°, 

and in addition, in the case of Neumann boundary conditions, assume that a = 1 and 6=0. Let Zt 
denote the number of zeros of /(•, t) in [0, 1]. Then 

(a) for t > , Zt is finite. 

(b) if (xq, to) is a multiple zero of f , then for all t\ < to < £2 we have Z^ > Zt 2 . 

By change of coordinates it is possible for an equation of type (6. 1) to be reduced to an equation in 
which the coefficients a = 1 and 6 = (see [3, 13]). 

The non-linear case has been studied by Galaktionov [13] and Angenent has obtained similar results 
in [3]. As the non-linear case is important to us we will state it briefly. Consider 

(6-2) ft = L(x,t,f,f x ,fxx) in Q, 

where L(x, t,p, q, r) is nondecreasing relative to the last argument r£l (the parabolicity condi- 
tion). Let & denote the family of all continuous solutions to (6.2) and let Bc^, We define the 
difference w(x,t) = f(x,t) — V(x,t), where / € & and V G B and the number of intersections 
of V as the number of sign changes of w(x, t) : 

Int(t,V) = Z t (w). 

If L is sufficiently smooth , the difference w satisfies a linear parabolic equation 

(6.3) w t = aw xx + bw x + cw , 
where the coefficients are given by Hadamard's formulae 

(6.4) a= [ L r (x, t, /, /*, 9f xx + (1 - 9)V xx )d9 > , 

Jo 

(6.5) b= [ L q (x,t,f,9f x + {l-9)V x ,V xx )d9, 

Jo 

(6.6) c= [ L p {x,t,9f + {l-9)V,V x ,V xx )d9. 

Jo 

If these coefficients satisfy the conditions of Theorem 6.1 then we can draw the same conclusions 
about the number of zeros of w. 

6.2. Sturmian Theorem in Mean Curvature Flow and Volume Flow. Using the Sturmian theo- 
rem Altschuler, Angenent and Giga [1] proved the following: 
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• The number of necks of an axially symmetric surface evolving by mean curvature is finite 
and nonincreasing for < t < T ([1] Lemma 4.7). 

• The necks converge to points on the axis of rotation as t — > T ([1] Lemma 5.1). 

• If the two limit points of nearby necks are distinct, then for any compact interval between 
the two limit points, there exists a uniform e > such that p(x±,t) > e in that compact 
interval for all t G (0, T) ([1] Lemma 5.2). 

• The limit \\vn t ^T p{x\,t) = p(xi,T) exists, and p{x\,t) converges uniformly to p(x±,T) 
ast -+ T and p(xi,t) G C°° (R x (0,T]) as long as p{x±,t) > ([1] Lemma 5.3). 

The first author ([5] Corollary 1.2, Lemmata 2.3, 2.4, 2.5 ) proved the equivalent for volume pre- 
serving mean curvature flow. 

6.3. Zeros of p" and ( — )• We use the Sturmian theorem to prove that the zeros of p" and of — 
are discrete and nonincreasing. Furthermore, we show that the limits of the zeros of p" and ^ exist 

as t -> T. 



Corollary 6.2. Assume {M t } to be a family of smooth hypersurfaces evolving by (1.2). Assume 
in addition that p(xi,t) > e,e > , for a < x\ < b ,t G (0, T"),T' < T. Then the set 
Z t (p") = < x\ G R : p"(xi,t) = g-§ = 0> is a discrete set in [a, b], for all t G (0, T'). Moreover, 
the number of zeros of p" is a nonincreasing function of time. 



Proof. We know from Lemma 4. 1 

d p" n — 1 , r -r 

-tip = Tiri2 + h \n- + P ■ 

at 1 + p' A p 

Differentiating this equation with respect to x\ twice we find that 7] = p" satisfies 

hr] 8rj 2 p' 2 hp' 2 f] \ 2(n-l)p' 2 

"r T-i I TTTT^ ~ T~T V 




VTT7 2 - (i + p /2 ) 3 (i + p' 2 )l) ' p 3 

As this is a nonlinear equation, we resort to the intersection comparison method discussed above. 
Suppose we have two axially symmetric surfaces evolving by mean curvature while preserving their 
respective volumes, then we can define the difference w(x\, t) = p'{{x\, t)— p2(xi, i) , where p±,p2 
are the respective radius functions of the first and the second surfaces. For pi > e > 0, i = 1, 2, 
we have bounds on the second fundamental form and all its covariant derivatives, thus we have 
smoothness. The difference w{x\,t) satisfies (6.3), where the coefficients are computed as in (6.4), 
(6.5), (6.6). As pi > e we have bounds for all the derivatives of pi from ([5] Lemma 2.5). Thus 
our integrands are bounded. The coefficients satisfy the conditions of the Sturmian theorem as the 
integrands are bounded, and thus we can apply Theorem 6. 1 and conclude that the intersections are 
discrete and their number nonincreasing in time. In particular if the function w has a multiple zero 
at (xo,to) then for all < t\ < to < ^2 < T' , the strict inequality Z tl (w) > Z t2 {w) holds, so 
the number of intersections of the two surfaces, Z t (w) is strictly decreasing at t = to- Finally as 



1 8 MARIA ATHANASSENAS AND SEVVANDI KANDANAARACHCHI 

pi = c is a solution of (6.7) where c is a constant, we can conclude that the zeros of p'[ := p" are 
discrete and nonincreasing in time. □ 

Corollary 6.3. Assume {M{\ to be a family of smooth hypersurfaces evolving by (1.2). Assume 
in addition that p(x\,t) > e, e > 0, for a < X\ < b ,t £ (0, T') ,T' < T. Then the set 

Z t ( (f) j = i xi e R : (f) (xi,t)=o\ is discrete in [a, b], for all t G (0, T'). Moreover, 
the number of zeros of ( — I is a nonincreasing function of time. 



Proof. From Lemma 4. 1 we obtain 



2 o / / rr\ 2 



d \\( H \ 2 /t7 t-,( h \\ 2 (P~k)H ( Hq 2 k 



dt J \P J p\ \P J / P V P 2 P 



p 



2 



As the Laplacian and the gradient in the above equation are computed on the hypersurface, we need 
the equivalent equation on R: 

(6.8) 

d ff H Y\ 1 d2 ff H Y\ ( P'p" , P' \ d ff H ^^ 



(n — D- 
dt\\pj ) (l + p' 2 )d 2 x i y\ P J J \(l + p' 2 ) 2 K 'p{l + p' 2 )J dxi \\p 



p (1 + p' 2 ) dx\ dx\ \\pJJ P V P J P 



p \pj 



.,2 



d fH^ x 2 



1 + p' 2 \dxi \ p 

As this is a nonlinear equation, we resort to the intersection comparison method discussed above, 
and as in Corollary 6.2 by comparing M t to a cylinder, where ( — ) is a constant and k = 0, we 



obtain the desired result. □ 

After discarding an initial section of the solution, without loss of generality we may assume that 
the hypersurface has I local minima of p', and I + 1 local maxima of p'. Let the local minima of p' 
be located at {t-j(t)}i<j<i , and the local maxima of p' be located at {Xj(t)}i<j<i+i, so that 

a < i-i(t) < i 2 {t) < ■■■ < n(t) < b, 

a<Xi(t) <X2(t) <■■ <Xi+i(t) <b. 

Similarly we may assume that the hypersurface has m minima of ( -^ J . Let the minima of ( -^ J 
be located at {£,j(t)}i<j< m , so that 

a<£i(t) <6(*) < ••• <Zm(t) <b. 

Lemma 6.4. (Convergence of zeros of p"). The limits 

lim ij (t) = ij (T) and lim X j (t) = Xj (T) 

exist. 

Proof. The proof is exactly the same as in ([1] Lemma 5.1). For the sake of completeness and as 
we use it in Lemma 6.5 we include it here. 
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Assume that ij(t) does not converge as t — > T. Then liminf t _>/r tj(t) < lim sup t _^ T ij (t) , and 
we can choose an xo £ (liminf f ^x 6 J '(t),limsup t _ >T ij(t)). Since ij(t) is continuous, there is an 
infinite sequence of times tk — > T at which ij{tk) = xq, and at which therefore p"(xo,tk) = 
holds. 

Consider the family of curves p'(t) on [a, b] and p'{t) obtained by reflecting p'(t) about the hyper- 
plane x\ = xo- Here p'(x\,t) = p'(2xo — x\ , t) , and p' is defined on [2xo — b , 2xq — a]. The 
curve p' corresponds to ^-p(2xQ — x\,t). From ([5] Corollary 1.2) we know that the zeros of 
p' for x\ € [a, b] and of pi for xi e [2xo — b , 2xo — a] are finite and nonincreasing. Therefore 
the zeros of w = p' — p', where x\ G [max (a, 2xo — 6) , min (6, 2xo — a)], are finite and non- 
increasing as well. Moreover, the number of zeros of w drops whenever p 1 — p' has a multiple 
zero. If p' — p' has a multiple zero at (xo,t*), U < T, that means if p'(xo,i*) = p'(xo,t*) and 
p"(xo,t*) = p"(xo,U), then p'(xo,t* + S) — p'(xo,t* + S) has at most one zero in the interval 
[xo — e , xo + e] (see Theorem 6.1). As p"(xo,tk) = —p"(xo,tk) due to the one being the reflection 
of the other and as p"(xo,tk) = 0, we have p"(xQ,tk) = making p"(xo,tk) = /5"(xo,ifc). As 
p'(xo,tk) = p'{xo, tk), this would mean that multiple zeros of w exist at (xo, tk) for tk > t*. This 
is a contradiction. We must therefore conclude that the tj (t) converges after all. 
The same argument also shows that Xj(t)'s converge. □ 

Lemma 6.5. (Convergence of zeros of -^). The limits 

lim^-(i)=^(T) 
exist. Moreover, the zeros of — converge ast—tT. 

Proof. The proof is exactly the same as in Lemma 6.4 after the following changes. Instead of 
the curves p' and p' we work with the curves ( — J defined on [a, b] and ( -j J (xi , t) : = 

— ) (2xo — xi , t) defined on [2xo — b , 2xo — a] . Instead of ([5] Corollary 1.2) we use Corollary 



6.3 above to say that the zeros of w = ( — I — I -S- I are finite and nonincreasing. The remainder 
of the proof is the same as in Lemma 6.4. Thus £j (t) converges as t — > T. As the zeros of — are a 

subset of the minima of ( — ) we conclude that the zeros of — converge as t — >■ T. □ 



6.4. Implication of the behaviour of the zeros. 

Lemma 6.6. Let H denote the mean curvature of the evolving hypersurface Mp If \VH\ < cin a 
closed region S C £1 x [0, T], then 4— < cin S as well. 

Proof. For the magnitude of VH we obtain 

\VH\ 2 = \V T1 H\ 2 + ■ ■ ■ + \V Tn H\ 2 . 

As M t is an axially symmetric surface, the mean curvature H is constant on the n — 1 dimensional 
sphere for a fixed x\ coordinate. Here we let x(Z, t) = x(xi, 9±, ■ ■ ■ , n _i, t), and for 2 < i < n 
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so that 
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IV-ffl 



\V T ,H\ 2 



d 



H 



As \VH\ is bounded in S, we have the same bound for I j^- 1 as well. 



i:lr i 



□ 



i = T 



t = 




Figure 4. The paths of # = d and # = C 2 in Q" 

We recall that fi" = {] t<T &" = Ut<r{ x ( / ' f ) : H {l,t) < c 2 - 5 : S > 0} . Let us define the 
following: 

Definition 6.7. In a connected component of Cl" , consider any two paths where H = C\ and 
H = C 2 such that < C\ < C 2 < c 2 - 5 , 5 > 0. (Recall that < c 2 < /i(t) < c 3 .) Le? 

Zi(t) = {/GM n :iJ(Z,i) =Ci}, W xi(Zi(t), t) = (x(li(t), t),h) , 

/ 2 (t) = {/GM n :F(/,t)=C 2 }, and Xx(l 2 (t),t) = (x(l 2 (t),t),h) , 

a(t) = min{xi(^i(t), t),Xi(/ 2 (t), t)} , and f3(t) = max{xi(7i(£), t),xi(l 2 (t), t)} . 

Here U(t),i = 1,2 is the curve in M n x [0, T) that parametrizes H = d and xi(U(t),t) the 
corresponding x\ coordinate. 

Lemma 6.8. With the above notation, there exists a constant c such that \xi(li(t),t)—xi(l 2 (t),t)\ > 
c> for all t < T. 

Proof. From Remark 5.2 we know that y |^;/> y |r-„ > e > 0. As the height is always positive in 
Q" we have p(x 1 ,t) £ C°°(R x [0, T\) in that region by ([5] Lemma 2.5). We note that this holds 
in Cl" even at t = T, as the height is strictly positive. Thus, there exists a constant C such that 
\VA\ \ft„< C for all t € [0,T] . As \VH\ 2 < n\VA\ 2 we have bounds for \VH\. From Lemma 
6.6 we know that there exists a constant d such that | J^ | < d in Vt" , for t £ [0, T]. Therefore 



Pit) dH 



«(*) 



dx\ 



dx\ 



< 



/ t. — dx\ < d ax i , 

Ja(t) UXi J a (fi 



(C 2 - d) 



\Htf(t),t) - H(a(t),t)\ < c'(P(t) - a(t)) , 

<\x 1 {h(t),t)-x 1 (l 2 {t),t)\, forall te[0,T}. 



□ 



We recall that the zeros of — are located at {£j(t)}i<j< m , so that a < £i(i) < £ 2 {t) < • • • < 
£m(t) < b , and that lim^y £j(t) = £j(T) exist from Lemma 6.5. 
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Lemma 6.9. For any fixed time t the one-dimensional Hausdorjf measure of £j (£) on the x\ axis 

satisfies J^ 1 (£j (t)) = ,for all i , ifj^f 1 (£j (0) ) =0 (i. e. £j (£) cannot develop any sudden thickness 
or horizontal length on the spacetime diagram [a, b] x (0, T] ). 

Proof. From Lemma 5.7, \A\ 2 is bounded in Cl' and, as p 2 < \A\ 2 , we have bounds for p in Cl'. 
From ([5] Lemma 2.5), the radius function p(x\,t) is smooth in the region Cl' t because of the lower 
height bound. J4? l (£,i{t)) > would imply a horizontal line on the spacetime diagram for a given 
time t. This would mean multiple zeros of — for the corresponding values of x\ on the axis of 
rotation. From Corollary 6.3 we know that the number of zeros of — is decreasing, therefore this 



cannot happen. 



□ 



In particular we rule out the scenario in figure 5. As the zeros of — can only develop in Cl' and as 
the height is always positive in this region, we have p(x\,t) G C°° (R x [0, T\) in Cl' . Thus even 
at t = T the zeros of ^ are well behaved. 



t = T 



t 




Figure 5. ^(^(T)) > 
For the next Lemma we recall that 

£t t = [x(l,t):H(l,t)>^}, and ft = \J U' t . 

t<T 
Lemma 6.10. There exists a constant c > , such that the one-dimensional Hausdorjf measure of 
the projection ofCl' onto the x\ axis satisfies Jff' 1 (I((l 1 )) > c for all t < T. 

Proof. To invoke Lemma 6.8, without loss of generality we choose 5 = |j|, C\ = ^ and 

Co 



^■. Hence Q,![ = (x(/,t) : H(l,t) < ^02}. From Lemma 6.8 we have \xi(h(t),t) 



xi(h(t),t)\ > c > , for all t < T. By the definition of (l' t and choosing k(t), i = 1,2, such 
that H(k(t),t) = Ct we know that x(k(t),t) g ftf t ,i = 1,2. As ^(/(jV)) > \xi(h(t),t) - 
xi (h (t) , t) I for any t, we have the desired result. □ 

Lemma 6.11. In any compact set ScO' where y > e > there exists a constant c(e) such that 
f > c(e) > 0. 

Proof. By definition we know that H > ^ in Cl' for all t < T. In any compact set S where 
y > e > , by Proposition 5.9 we know that \A\ 2 is bounded and hence we get an upper bound for 
p. Therefore we have a lower bound for — in S. □ 

Therefore in Cl' , — can only tend to zero when y tends zero. We will investigate this in the next 
section. 
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7. The Singularity 

We break up the investigation of the singularity into different cases, depending on the value of |fe|. 

With the aid of an additional assumption which involves a lower height bound, we will prove that 

the singularity is of type I. From now on all our calculations are done in M 3 for two dimensional 

surfaces. 

First we will state and prove two useful results in a restricted setting. We consider a two dimensional 

surface in ]R 3 and follow a method as used by Huisken in [17]. 

7.1. The region S. Let St C Cl' t and S = [j t<T St- For this region S we assume that there exist 
constants c\i , C13 > such that 

\A\ 2 



< C12 and y\ r > c 13 . 



s 
S 



The following Lemma corresponds to Lemma 5.2 in [17]. 



Lemma 7.1. Under the above assumptions, there exists a constant C15 such that ^ < C15 in St 



Proof. From Lemma 4. 1 we compute the evolution equation for -^ 

d_(±\ A / <? \ , 2 /nrr w / <? \ \ , 9 ( I Jl _2 o,..^ , h 

dt \H, 

We know in S t 

|2 ^- „ tt2 o ;„„ ** 7,2 , 2 i \ |2 ^- „ rj2 



A (l) + ^( Vi/ ' v (/l)) + l (p 2 -9 2 -2M + |(fc 2 -M • 



i^r < ci2^ , -2&p < r +p = i^r < cia^ , 

ct h 2c-\ 

0<c 2 </i<c 3 , H > — , such that — < — . 

2 H c 2 



When £ > 2 J^f 2 we obtain 

p 2 -q 2 - 2kp + 4( fc2 " M < C 12^ 2 " 8-ci2# 2 + ci 2 # 2 + — (c 12 H 2 + ^ci 2 # 2 
il c 2 c 2 V 2 

= c 12 ^( 2 -^)<0, 



as ^ > 1 . Therefore, when -ft > 2-. / 2c ^ Cl2 , from the parabolic maximum principle and viewing 
the surface as periodic 



, , Q Q /2c 3 ci2 

— < max max — , max — , l\ 



H ~ \ So H 1 ds t H 1 V c 2 

t<T 



We recall that q = (v,ii)y 1 . As max&s t -ft < maxg^, -^ — and as y 1 \ds t ^ ~ > tne right 

t<T i<T C2 «T C13 

hand side in the above estimate is bounded. 
Similarly when jj < -2J 2c:i c Cl2 we have 

dt \HJ ~ \HJ H\ ' \H 



SINGULARITIES OF AXIALLY SYMMETRIC VOLUME PRESERVING MEAN CURVATURE FLOW 23 

Thus from the maximum principle we have ^ < cu . As - < c\ we obtain 

\q\ < c u H = cu{p + k) < a 5 p , 
as desired. D 

Remark 7.2. (i) Note that ^ is a geometric quantity that corresponds to the slope \p'\ of the gen- 
erating curve p. Therefore, Lemma 7. 1 gives us a gradient bound in the region S. 

(ii) Assuming that the singularity develops in S, there is a point on the generating curve that ap- 
proaches the axis of rotation as t — > T. By definition, as y\r s > C13, we have y\ds t ^ c v&f or all 
t < T. As the boundary of the domain has a lower height bound, and as the above gradient bound 
holds in that domain, St cannot collapse to a point as t goes to T. 

The following Proposition corresponds to ([17], Proposition 5.3). 

Proposition 7.3. If the singularity develops in the region S, then there exists a constant C > such 
that the second fundamental form satisfies 

,9 1 

max \Ar < C , 

S t ' ' - T-t 

for all t < T. 

Proof. From (4.3) we know 



^-y' 1 = (H- h)py- 1 >(H- c^py' 1 . 
at 

Asp 2 < |^4| 2 < cuH 2 , and using Lemma 7.1 

y- 2 = p 2 + q 2 < (1 + C 2 5 )p 2 , 
d -1 ( 1 \ 1 -2 



dt a ~ Vv 7 ^ ) 7(TT4)' _ 

As S contains the singularity, for t near T, we have p — > 00. For p > lyjcyicz 

1 1 

p - c 3 > - p . 



y/CY2 2y/c~Y2 

Therefore 

d -1 -> I 1 -2 ^ 1 - 

-j) > n r— v—F= .> = y > „ , — ,, , ttv 



df ~ 2^i 1 \/(l + c{ b ) ~2y/c^(l + c 



2 
15 J 



Let U(t) = max5 t y L . By renaming the constant — , 1 3-r = e we obtain 

jU(t)>eU*(t)^jU- 2 (t)<-2e. 

Since U~ 2 {t) tends to zero as t — > T, we integrate from t to T and obtain 

1 



U(t) = maxy < 



St V 2e ( T ~ t) 

As \A\ 2 < C12H 2 and H = k + p < c±p + p < (c\ + l)y _1 we get the result. □ 

These results will be useful when we consider the different cases outlined below. 
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7.2. Different cases. We consider different scenarios, depending on the value of \k\. For this we 
divide the hypersurface into different regions. We now state an additional assumption, which will 
only be used when | k\ is unbounded; it presumes a height bound on the boundary of Q\ independent 
of time. 

Assumption 7.4. We assume that there exists a cqq > 2 such that y|g^b > cfor some c > Ofor all 
t <T. 

The different cases we look at are shown in the figure below. 



The evolving hypersurface 



|fc| is bounded |fc| is unbounded 

on M t for all t < T Assume 7.4. 

Singularity develops in 



* 



n' n o b o' n n« 



Due to the lower height bound in Cl' ( Lemma 5.1), and Proposition 5.6, we know that \A\ 2 is 
bounded in Cl'. Therefore the singularity can only develop in Cl' . As M t = Q t U £l t , we have 
Cl' = ( Cl' n & ) |J ( Cl' n Jl" ) . Therefore the cases shown in the diagram below present all the 
possible conditions under which a singularity can develop. The case where |fc| is bounded on M t 
for all t < T is considered in section 7.3. The case where \k\ is unbounded and the singularity 
develops in Cl' n ST is considered in section 7.4 and the case where \k\ is unbounded and the 
singularity develops in Cl' n £$ is considered in section 7.5. The latter two cases are investigated 
under Assumption 7.4. 

7.3. Case I : | k\ is bounded on M t for all t < T. In this case it can be shown that '-jX- is bounded 
in Cl' as follows. As 

\k\ < Ci 6 , 
we get 

p 2 - 2c w p <k 2 +p 2 - 2c 16 p <k 2 +p 2 + 2kp . 

Hence 

|^| 2 k 2 +p 2 c 2 w +p 2 



< 



H 2 k 2 + p 2 + 2kp p 2 — 2ciqp 
In the case of a singularity developing, we can assume without loss of generality that p > 4ci6, 
such that p — 2ci6 > f ■ Therefore, for p > Ac\§, 

\A\ 2 k 2 +p 2 ^ c 2 % +p 2 ^^+P 2 _17 



< , U r < 



H 2 k 2 +p 2 + 2kp ~ p{p - 2ci 6 ) ~ vL 8 
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From Lemma 5.1 we know that infrv, y > d . By letting St ■= Cl' t , from Lemma 7.1 and Proposi- 
tion 7.3 we conclude that the singularity is of type I. 

7.4. Case II : \k\ is unbounded and the singularity develops in Cl' n &". Cases II and III are 

studied under Assumption 7.4. Let 

(Ci' t n nty = |x(i, t) e (Ci' t n n b t ) ■. k(i, t) > 0} and 
r n' n n^ =\J (Ci' t n njf 



Let 



(Ci' t n n£) = jx(z, t) e (Ci' t n n b t ) ■. k(i, t) < 0} and 
'Cl' n ft b ) = J (n{ n of 



In (0'n^) + , we have ^ = ^^^ < 1 as k > 0. 

As coo > 2 we have 2(e ^ 1} < 1 . Let S := 1 - J 2 ^-i) > °- In (^' n ^ b ) we have 1*1 = ~ k 

as k < 0. Therefore in ( fV n f^ ) , by definition of the domain ( see 4.3.4, ) 

--<i-s, 

P 
5p< H . 

This estimate combined with 

\A\ 2 = k 2 + P 2 <(^^ + l) P > = p^ P \ 
\2(coo-l) / 2(coo-l) 

gives us a bound for '-^r in ( (V n fiM .As '-jp- < 1 in ( Cl' n & J we have an upper bound for 

LL. in Cl' n f^. From Lemma 5.1 we know that infr-, y > d > and from the Assumption 7.4 
we know that y\ dn b > c > for all t < T. By letting St := Cl' t n r^ , from Lemma 7.1 we get a 
gradient bound in Cl' n r^ . By Proposition 7.3 we conclude the singularity is of type I in this case 
as well. 



7.5. Case III : \k\ is unbounded and the singularity develops in Cl' n fi". First we prove a 
couple of useful results in the domain Cl' HflK We show that the gradient function is bounded in 
Vl\ in the following Lemma. 

Lemma 7.5. There exists a constant c\i such that v < en in 0,1 for all t < T. 

Proof. From Corollary 4.3 we know 

^--a\v = -\A\ 2 v+p 2 v 3 - 2vk 2 {v 2 - 1) . 
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For v > a/coo we have v 2 — 1 > ( Coo ~ 1 ) v 2 , giving us 



('0(1 



d At„, s- l/1|2„, , „2„,3 o f C °° 1 1 l.2„,3 



— - A )v < -\A\ 2 v+p 2 v 6 - 2 ^ k 2 v 

at J V c oo / 

^OO- 1 ^ v 3 2 f _ COO A 1 ' 

. coo y \2( 

as — > . / „, Coo _ 1 \ in S7^- By the parabolic maximum principle 



Cqo / V 2 ( c oo- 1) P 2 



1*1 

p "" \l 2(coo-l) 

v < max max v , max v , v c 00 

\ *<T / 

From Assumption 7.4 we know that y > c > on df^ = dOJL From Lemma 5.3 we know that 
vy < C4 on M t for t < T. Therefore v < — on dVt\ for all t < T. So we have 

< ( C4 ^ 

f < max max!) , — , v c oo =: c\i ■ 

V < c J 



As infp., y > c' > and J/LqU > c > for all t < T (Assumption 7.4), and as the gradient is 
bounded in the domain Vi\ for all t < T , by an argument similar to that in Remark 7.2(h) we see 
that the domain Cl' n £l" does not collapse as t tends to T. 

Note that for this case (singularity develops in Cl' n fi" ) we do not have an upper bound for '-jfi- and 
thus we cannot use Proposition 7.3. Our approach here is to investigate the behaviour of the curva- 
ture along paths going to the singularity, i.e. smooth curves in the domain SI x [0, T) approaching 
the singularity. First we consider paths along which — stays bounded away from zero. 

Proposition 7.6. Let A t be a path going to the singularity in fl» n O'. If — > ci8 > on At for 
some constant c±s > 0, then there exists a constant C such that 

.,.21 C 

\A\ L < . 

1 I I At - T _ t 

Proof. The proof is along the same lines as that of Proposition 7.3 with only a few changes. From 
(4.3) we know 

-jtV" 1 = (H - h)py~ l > (asp - c 3 )py~ l . 
at 

As p — > oo on A t , for p > 2-^- we have c\gp — C3 > \c\$p ■ 
Therefore 

d -1 ^ c i8 2 -1 c i8 y~ 3 1 

—y > — p y = 7T- , as p = — . 

dt y ~ 2 1 y 2 v 2 ' vy 

As v < C17 in £l\ 

d _! Cis -3 

dt y ~ 2c 



-17 
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Let U(t) = y~ l I . .By renaming the constant ^§- = e we obtain 

I l\t ' " 17 

Since U~ 2 {t) tends to zero as t — > T on A t , we integrate from t to T and obtain 

1 



U(t)=y' 1 \. < 



VHT - 1) 



From Proposition 5.5 we know that '—± < max(ci , n — 1) in Cl' D Q%. As we are considering a two 
dimensional hypersurface, we have \A\ 2 = k 2 + p 2 < c\p 2 + p 2 < {c\ + l)y~ 2 , along any path 
A t in fV n ri" where — > ci§ > , and thus we get the result. □ 

In what follows we prove that a singularity that develops in $Y n £$ is of type I by way of contra- 
diction. We start by considering a statement which is the negation of a singularity of type I. Then 
we prove that this statement is false by using a rescaling procedure similar to that used in [18]. 

Statement 7.7. For all c > there exists t(c) € [0, T) , such that 

max \A\ 2 (l,t(c)) > — . 

Proposition 7.8. Statement 7. 7 is false in Q! n fi". 

Proof. Suppose Statement 7.7 were true in Cl' n O". Let us choose a sequence {q}j G n such that 
linij^oo Ci = oo. If there exists more than one £(cj) for a particular c» , we choose the earliest 
t(ci) =: ti such that the above inequality is satisfied. Then 

mac\A\ 2 (l,ti)(T-ti) >Ci. 

The sequence {a} going to infinity forces the left hand side of the above inequality to also go to 
infinity. Therefore maxj e M n |^4| 2 (^, U) has to go to infinity as a goes to infinity. That means that 
ti approach T as i goes to infinity. 
We now rescale Cl' n f^ as follows: Let 

ctj = max \A\(l, t) 
t<u 

and the points k be defined by the equation 



\A\(li,U) = max\A\(l,U) . 



Then we have 

(7.1) Oi>\A\{li,U)> 



T-U 
We consider the family of rescaled surfaces ^#j iT defined by the following immersions: 

(7.2) Xj(-,r) = Oj (x(-,a 4 " 2 r + £j) - (x(/j,tj),ii)ii) , 

where r € [— afti, 0] . For this rescaling we have 

H i (;T) = ar 1 H(-,ar 2 T + t i ) and \A\\(; t) = aT^A^^a^T + U) . 
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1*1 

Cl'nQK Therefore 



This rescaling guarantees that \Ai\ < 1. From Proposition 5.5 we know that '—± < max(ci, 1) in 



\A\ = y/k 2 + p 2 < c 5 p < c 5 y x , 
where C5 = \J\ + (max(ci, l)) 2 . Hence we obtain 

(7-3) \Ai\ = a^ 1 \A\ < a^ l c 5 y~ l = c 5 (ajy)~ 1 = c 5 y~ l . 

Thus the rescaled surfaces do not float away to infinity. As | A4 \ < 1 for all i, by employing a diago- 
nal sequence argument, we find that the rescaled surfaces ^#^0 converge to a smooth hypersurface, 
which we denote by ^#0 . We explain this process in the next paragraph. 

Without loss of generality let us assume that the singularity develops at the origin. We translate 
each generating curve to have (x(7j, ti) , i\) = 0. We rename them, using the notation for the radius 
function, p, and work with them for the rest of the section. As Ai(li, 0) = 1 , from (7.3) we have 
y(h, 0) < C5 for all i. Thus we find a converging subsequence of points. We look at the unit tangent 
vectors at x\ = of the corresponding subsequence of curves. By identifying these unit vectors 
with points on the sphere, we observe that a subsequence of these points converges and so obtain 
a subsequence of tangent vectors that converges. As A\ < 1 on the hypersurface, we can roll a 
unit radius ball anywhere on the hypersurface. Thus, after rotating each curve of the earlier subse- 
quence such that all the unit tangent vectors at x\ = are horizontal, we have a uniform gradient 
bound within a small radius \x±\ < S < 1. Therefore by Arzela-Ascoli we find a subsequence that 
converges uniformly within the interval \xi\ < 5, in which the uniform gradient bound holds. Re- 
peating this process along the curve we can find a uniformly converging solution within any given 
compact set. This process is explained in detail in [20] Chapter 5. 

Now we consider paths that lie in &' n fi" in the spacetime diagram £} x [0, T) . We look at the value 

of — on these paths. We note that as - < c\ on M t for all t < T (Proposition 5.4), — < 1 + c\ on 

M t as well. Thus f is bounded on M t for all t < T. 

From Corollary 6.3 we know that the zeros of — are discrete and non-increasing with time. And 

from Lemma 6.5 we also know that the zeros of — converge as t tends to T. Thus the zeros of — 

are well behaved. 

The quantity — can be zero either when H = or when p = 00. However, the zeros of — that 

correspond to the zeros of H lie outside Cl'. In conclusion the only point where — can tend to 

zero in the domain Cl' n fi" is the singular point. That means, if we follow any path in the domain 

Cl' n O", the value of -^ cannot be zero as long as we are not at the singularity. 

We consider at different paths in Cl' n fF for t G [0, T). We categorize these paths by the behaviour 



H_ 
P' 

paths. 



of — , and show that the sequence Hi converges to zero as i goes to infinity, along each of these 



(1) Paths going to the singularity - paths where ^ > c > for some positive constant c 

From Proposition 7.6 we know that on these paths \A\ 2 < 7^-^ for some positive constant C. 
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Therefore on these paths 



H 2 <n\A\ 2 < " C 



From (7.1) 



As -?^ < 1 for t < tj , 



H 2 = a~ 2 H 2 < 



T-t 
T-tA ( nC 



Ci J \T — t 



nC 
H 2 < , for t < U 



Ci 

As Ci goes to infinity, Hi converges to zero on these paths. 



(2) Paths going to the singularity - paths where 3 _ >. o as t — > T 

As a.i > p for t < tj, on these paths we have 

Hi = a7 1 H < >0. 



On these paths Hi converges to zero as well. 



(3) Paths not going to the singularity 

On these paths \A\ 2 < c for some constant c, that depends on the end point of the path. As c^ goes 
to infinity, we have 

Hi = aT x H < a r l ^\A\ < a^y/nc -)• . 
Therefore on all paths in Cl' D fi" , Hi converges to zero as i goes to infinity. The limiting solution 
^#o is a catenoid, as it is the only axially symmetric minimal surface with zero mean curvature. 

We are now in a position to show that we have a contradiction: In order to get a better understanding 
of the original surface we rescale back ^i$ for large i, and show that the estimate vy < C4 would 
not hold on that (the original) surface. 

We denote the quantities associated to the catenoid ^#0 by a hat ". We obtain the catenoid ^#0 
by rotating y = C5 cosh(c^" 1 xi) around the x\ axis, where x 1 is the x\ coordinate of the limiting 
surface ^#o- F° r any ei > and for any l§ G M n we have 

\v{lo)y{h) - Vi(l , 0)yi(lo, 0)| < ei for large i . 



For the catenoid v = yl + y' 2 = Wl + sinh 2 (c 5 x\) = cosh(c 5 x\). As yi = aiy and Vi = v 
we have 

c 5 cosh 2 (cg" 1 xi(/ )) - ei < aiv(l ,ai 2 T + ti)y(l ,a^ 2 T + ti), 
(7.4) ^- (cosh(2 C 5 1 x l ) + l) - — < vy for % > I . 

£Ct% Oii 
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As 

xi{l ) = (x(/ ),ii) = lim aij ( Xi(/ , oJ 2 t + tj) - sy ) 

where xy := (xj, i^ , for any e 2 > we have 



Xi(/ ) - aj (xi(Iq , a^ 2 r + tj) - Xij) 



< e 2 for large j , 



(7.5) aj [xi(l , a- 2 r + tj) - x Xj j - e 2 < xi(l ) < e 2 + aj (xi(fo , «j 2 r + tj) - xyj , 
for j > Jo- For given ei, e 2 and for a fixed i > Iq , we pick large j > Jo , j » i, such that 

«j v V c 5 V a iJ J J 

with e sufficiently small and with C4 as in Lemma 5.3. Then we choose Iq G M n such that 

(7.6) xi [l , aj 2 T~o + tj) > xij + I 

holds. Considering M t as a periodic surface, we can find points on the hypersurface that lie an e 
distance away from x\j. Therefore 

1 / /4aj / t\ 



(7.7) xi U 0l a j T + tj I > xij H I e 2 + 2c 5 log I I c 4 H 1 

We will use the inequality (7.7) to obtain a contradiction. By (7.5) and (7.7) 

Xi(l ) > aj \xi(l , aj 2 T + tj) - xijj - e 2 > 2c 5 log ( — - ( c 4 + — j - 1 j > . 



Therefore 



_£5_ 
2a., 



cosh I — (aj(xi - xij) - e 2 ) j < — - cosh ( —xi(l ) j for j > J . 



By (7.4) 

S ( C ° Sh (2^ {a ^ Xl ~ X ^ ~ €2) ) + "I < 5" ( COSh (d^* ) + " a7 < ^ ' 

for ei,e 2 ,i,j and /o as previously chosen. As aj (x\(lo,aJ 2 T + tj) — xij) — e 2 > 0, by the 
choice of lo 

(7.8) -^_ f e 2^K(^-^i)-^) + l) _ £1 < W y . 

4a« V / «j 

But from (7.7) we have 

(7.9) _£L |' e 2^( Q! j( a: i(*o,a7 2 -n3+*j)-xi J -)-e2) + 1 N \ _ £1 > C4 _ 

4a« V / ai 

From Lemma 5.3 we know that wy < c 4 . Therefore (7.8) and (7.9) contradict Lemma 5.3: by ex- 
amining the rescaled surfaces we find that the estimate vy < c 4 does not hold on the corresponding, 
non-rescaled, hypersurfaces near the singular time T. 
Therefore we have a contradiction and Statement 7.7 is false. □ 
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Hence there exists a constant c > such that for all t G [0, T) , 

max \A\ 2 (l,t) < — ^— . 

leMn ' ' v ' ' ~ T-t 

We have proved that the singularity is of type I if it develops in fi' n fi". The combination of cases 
I, II and III, gives the proof of Theorem 1.1. 
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